Experimental Realization of Optimal Noise Estimation for a General Pauli Channel 
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where i,j= 1, ..., d, & and j denote the estimated and the true 
values for the 7-th component of the noise parameter respec- 
tively, Eg denotes the expectation with respect to the measure- 
ment procedure employed, and M is a projective measurement 
operator. 

In this work we present an experimental implementation of 
an optimal quantum channel estimation scheme for a qubit 
Pauli channel (PC). The action of such a family of channels 
on the density operator p of a qubit can be described as |4|] 
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We present the experimental realization of the optimal estimation protocol for a Pauli noisy channel. The 
method is based on the generation of 2-qubit Bell states and the introduction of quantum noise in a controlled 
way on one of the state subsystems. The efficiency of the optimal estimation, achieved by a Bell measurement, 
is shown to outperform quantum process tomography. 

PACS numbers: 42.50.Dv,03.67.Bg,42.50.Ex 

Introduction.- Quantum noise is unavoidably present in 
any realistic implementation of quantum tasks, ranging from 
quantum communication protocols [1] to quantum informa- 
tion processing devices and quantum metrology (i ISfl . The 
performance and the optimization of quantum tasks quite of- 
ten depend on the level of noise which is present in the phys- 
ical realization considered. It is therefore of great interest to 
develop experimental methods to estimate the level of noise in 
the system under examination as precisely as possible. Quan- 
tum process tomography (QPT) [4], which has already been 
implemented in various experimental realizations |5l|g]» rep- 
resents a well-known method to identify an unknown noise, 
but it lacks the notion of efficiency. In many realistic scenar- 
ios, however, some a priori information on the kind of noise is 
available and therefore the problem of measuring it is equiv- 
alent to estimate few noise parameters in the most efficient 
way. This is the context of quantum channel estimation |7|], 
which is based on quantum state estimation theory, a merge of 
statistical estimation theory and quantum physics pioneered 
by Helstrom and Holevo ISl |_10G • The aim of this Letter is to 
provide the first genuine experimental application of quantum 
channel estimation theory. The experimental realization pre- 
sented here is based on a quantum optical setup, but it opens 
new perspectives of applications to a great variety of physi- 
cal scenarios and quantum technologies, from atomic to solid 
state systems. 

A quantum channel estimation problem is generally formu- 
lated as follows. We need to estimate the true value of the d- 
dimensional parameter for a given smooth parametric fam- 
ily Yq of noisy quantum channels. The estimation scheme is 
twofold: first we prepare a quantum system described by the 
density operator p as input to the channel Tg, then we perform 
some quantum state measurement on the output state Tg (p) 
in order to estimate the channel parameters in the most effi- 
cient way. Thus, the problem is to seek an optimal input p for 
the channel and an optimal measurement on the output state 
Tg (p). The notion of optimality is here based on the mini- 
mization of the covariance matrix Vg[p, M] 
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where cr is the identity operator, {<x;} (i = 1,2,3) are the 
three Pauli operators <T x ,o- y , cr z respectively, and {/?;} repre- 
sent the corresponding probabilities (2?=o Pi = !)• The family 
of Pauli channels represents a wide class of noise processes, 
that includes several physically relevant cases such as the de- 
polarising channel, which will be considered in the following, 
the dephasing and the bit-flip channels. 

The optimal channel estimation scheme is achieved as fol- 
lows [9fl. The optimal input state is represented by a Bell 
state for two qubits, for example the singlet state \\ff~) = 
4=(|01) - 1 10)), where only one of the qubits is affected by the 
noisy channel while the other one is left untouched. The opti- 
mal measurement consists of a Bell measurement on the two 
qubits at the channel output, namely the projective measure- 
ment M = {|r><rUr><^ + U0"X0"U0 + ><0 + l}. The out- 
come probabilities then provide an optimal estimation of the 
channel parameters p\. 

As mentioned above, this scheme is optimized by minimiz- 
ing the covariance matrix of the estimation error (Q]). Accord- 
ing to the quantum Cramer-Rao theorem JlOfl , the minimum 
covariance matrix in this case is given by |9|] 
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where J p [pme] is the quantum Fisher information matrix | 
of a maximally entangled input state pme- We want to point 
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FIG. 1. (a) Experimental setup. Photons A and B are spatially and 
temporally superimposed on a symmetric beam splitter (BS). The 
optical path delay Ax allows us to vary the arrival time of the pho- 
tons on the BS. Photons are collected by using an integrated system, 
composed by a GRIN lens (GL) and a single mode fiber, and then 
detected by single photon counters. The same setup allows us to per- 
form the ancillary assisted quantum process tomography (AAQPT) 
after removing the BS. Quantum state tomography (QST) 1 19] on 
the output state is performed by using quarter- wave plates (QWPs), 
half- wave plates (HWPs), and polarizing beam splitters (PBSs). (b) 
Scheme of the implemented Pauli channel. t\, h,h represent the time 
intervals of cr x , cr y or cr z activation. Both t\, t 2 , h and the repetition 
time T can be varied by a remote control. 



out that this scheme is optimal for any number of input qubits. 
Actually, no additional entanglement among the input qubits 
and no collective measurements at the output can increase the 
efficiency of the present scheme [91] . Moreover, it can be also 
straightforwardly generalized to estimate any general noise 
process of the form 0, where the cr operators are replaced 
by any set of unitary operators V* such that Tr[VjVt] = 2Stj. 
The same scheme can also straightforwardly extended to es- 
timate any generalized Pauli channel for quantum systems in 
arbitrary finite dimension Q2D. 

We will now present the experimental implementation of 
this optimal estimation scheme for a quantum optical setup, 
where the state of the two qubits is represented by polariza- 



tion states of two photons and the action of the Pauli channel 
is introduced in a controlled way by employing liquid crys- 
tal retarders, as explained in the following. The method has 
been first applied to estimate a general Pauli channel, with 
independent values of the probabilities pt. Then it has been 
applied to a depolarizing channel (DC), namely the case of 
isotropic noise, with p\ = pi = P3 = f , where the parameter 
p completely specifies the channel itself, and the minimum 
variance (Eq. (Q} for the one-dimensional case) is given by 
p(l - p). In this case the procedure simplifies and, in the 
following, we show that only two projective measurements, 
M' = {\i//-)(i//-\, 1 - |<A"X*A"|}, are needed. 
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FIG. 2. Experimental probabilities of measuring the four Bell states 
obtained for four different cases of anisotropic noise. The black 
boxes report the corresponding theoretical values, a) Identity: noise- 
less channel, b) <r y : only one Pauli matrix, <r y is acting on the state 
\i//~). c) Partially anisotropic DC: cr x and cr z operate for the same 
time interval, in fact the probabilities of measuring the states \if/ + ) 
and |0 + > are equal, d) Totally anisotropic DC: each Pauli operator 
operates for a different time interval. 

Experimental scheme. - Different techniques have been ex- 
ploited to experimentally implement a PC acting on a sin- 
gle qubit state Ill2l4l5ll . The optimal noise estimation pro- 
tocol, proposed in this work, was implemented by the inter- 
ferometric scheme shown in Fig[T^). Precisely, a two-photon 
entangled source 111 ill generates the two-qubit singlet state 
\\ff~) = -j=(\HV)ab - \VH)ab), where two qubits are encoded 
in the polarization degree of freedom, with H (V) referring to 
the horizontal (vertical) polarization of photons A and B. In 
our setup, the single qubit noisy channel is operating only on 
one of the two entangled particles (i.e. photon A). The general 
Pauli channel (PC) consists of a sequence of liquid crystal re- 
tarders (LCI and LC2) in the path of photon A. The LCs act as 
phase retarders, with the relative phase between the ordinary 
and extraordinary radiation components depending on the ap- 
plied voltage V. Precisely, V n and V^ (FigHh) correspond to 
the case of LCs operating as half-wave plate (HWP) and as 



a) Optimal Protocol 



Pexp 




0,0 0,2 0,4 0,6 0,8 1,0 



5/T 



p 



1,0 
0,8 
0,6 
0,4 
0,2 



0,0 



b) AAQPT 
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FIG. 3. Noise parameter estimation for the DC case, (a) Measured values of p exp vs | by implementing the M' projective measurements. 
Continuous red line corresponds to the theoretical behavior. Inset: experimental values of the standard deviations for the optimal protocol 
implemented by the M' projective measurements. They are obtained by propagating the Poissonian uncertainties. The solid line represents 
the expected theoretical behavior, b): Experimental probabilities associated to the experimental matrix^ vs |. Values of p are obtained by 
maximizing the fidelity F between theoretical and experimental matrix^. Error bars are calculated by considering the Poissonian uncertainty 
associated to the coincidence counts, and simulating different matrices of the process, obtaining, in this way, different values of p. Inset: 
experimental values of the standard deviations for the AAQPT. These have been calculated by a simulation based on the poissonian uncertainty 
associated to the coincidence counts. The solid line represents the optimal bound. (c),(d) Experimental (left side) and theoretical (right side) 
matrices \ f° r T - 6 and T & 36. 



the identity operator, respectively. The LCI and LC2 optical 
axes are set at 0° and 45° with respect to the V-polarization. 
Then, when the voltage V„ is applied, the LCI (LC2) acts as a 
o~z (°~x) on the single qubit. We were able to switch between 
V^ and V n in a controlled way and independently for both 
LCI and LC2. The simultaneous application of V n on both 
LCI and LC2 corresponds to the cr y operation. We could also 
adjust the temporal delay between the intervals in which the 
V„ voltage is applied to the two retarders. We define t\, t 2 , h 
respectively as the activation time of the operators <j x , cr y or 
cr z and T is the period of the LCs activation cycle, as shown 
in Fig [it). 

Experimental implementation of the Pauli channel 
(anisotropic noise). - A general PC was generated by varying 
the four time intervals t\, t 2 , h and T. The intervals t t are 
related to the probabilities p t (i = 1, 2, 3), introduced in Eq. 
©, by the following expression: p t = j. The probability 
Po of the identity operator is given by po = 1 - | (with 
S = t 1 + t 2 + h). 

To obtain the probabilities associated to the four projection 
operators M, we measured the coincidence counts between the 
two outputs of the BS. In fact, these probabilities are related 
to the interference visibility measured by the interferometer in 
FigHt). The half- wave plate (HWP*) and quarter- waveplate 
(QWP*) of Fig Ok) were used to project the noisy state onto 
the four different Bell states. 

Different configurations of the noisy channel were investi- 
gated by implementing the optimal noise protocol estimation 
for each configuration. A summary of four relevant experi- 
mental results, corresponding to different probabilities associ- 
ated to the Bell states, are given in Fig El 

In the measurements shown in Fig El case (a) correspond to 
a noiseless channel (identity transformation) while the cases 



b), c) and d), correspond to different complete noisy channels 
with po = (i.e. we set T = S). For each process, the first 
column shows the relative weights between the Pauli opera- 
tors acting in the channel. From these values it is possible to 
calculate the theoretical ones. For instance, let us consider the 
process d) where the cr z , cr y and cr x act respectively for |, ^ 
and ^. The expected values of/?/ are, for this process, po = 0, 
Pi = \*Pi - | and/?3 = |. The slight disagreement between 
the expected theoretical values and the experimental measured 
ones are mainly due to the finite rise and decay times of the 
electrical signal driving the LC devices. 

We have implemented the protocol by using always the 
same input state and projecting it on the Bell basis. It is worth 
noting that this is totally equivalent to entering the PC with 
the four Bell states and to projecting them into the \if/~) state. 

Experimental implementation of the depolarizing channel 
(isotropic noise). - The condition t\ - t 2 = h corresponds to 
the depolarizing channel, with the three Pauli operators acting 
on the single qubit with the same probability p - | = h+t ^ +h . 
This parameter was changed by fixing the times t; and varying 
the period T. The optimal protocol to estimate the value of p 
was realized by using the Bell state \if/~), as mentioned above. 

The DC was activated on photon A. In this case the pro- 
jective measurement M' - {\if/~)(if/~\, t - \if/~)(if/~\}, consist- 
ing of just two projectors, is sufficient to optimally estimate p 
and has been performed for several noise degrees. For each 
level of noise, we estimated the channel parameter p exp as 
Pexp = N + c . where d nt are the coincidences between the 
two outputs of the BS in interference condition and N ss is the 
number of events in which the two photons are detected on 
the same BS output side. N ss was estimated by knowing the 
amount of coincidences out of interference. The typical peak 
interference measured for the state 1^") as a function of the 
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path delay Ax is shown in Fig H of the Supplemental Material. 
In Figl2k) we report the experimental values p exp correspond- 
ing to the different values of T. In the corresponding inset 
we show the p exp errors evaluated by propagating the C[ nt and 
N ss Poissonian errors. They are in good agreement with the 
expected theoretical behavior. 

Ancillary assisted quantum process tomography. - The ex- 
perimental results, just discussed for the optimal estimation 
of the depolarizing channel, have been compared with the 
probability values of p which can be obtained by exploiting 
the Ancillary assisted quantum process tomography (AAQPT) 
Ill4 Il6l4l8ll . The action of a generic channel operating on a 
single qubit can be written as 8[p] = Tth=oXij°~iP°~ji where 
the matrix Xij characterizes completely the process. 

AAQPT is based on the following procedure: i) prepare 
a two-qubit maximally entangled state and reconstruct it by 
Quantum state tomography (QST) Ill9ll ; ii) send one of the 
two entangled qubits through the channel 8; iii) reconstruct 
the output two-qubit state by QST and obtain, in this way, the 
matrix Xij from the two-qubit output density matrix. For a DC, 
the matrix xij is expressed as Jj] 



(4) 



We implemented the AAQPT algorithm by injecting the 
state \ifr~) into the DC and we reconstructed by QST the den- 
sity matrices of the input and output states for several noise 
degrees [see FigHk)]. We obtained the experimental matrix 
Xexp for different values of T and, for each value of T, we 
found the parameter p maximizing the fidelity between the 
experimental Xexp and the theoretical^^ process matrices. 
The experimental results are shown in Fig Ob). Even in this 
case the theoretical behavior is fully satisfied. However, com- 
paring these results with those obtained by the optimal pro- 
tocol, we observe that the latter leads to the same results, but 
with a much lower number of measurements. In fact, in this 
case, only the two projections M' are needed while, to imple- 
ment the AAQPT algorithm, 16 measurements are necessary. 
Moreover, by adopting our experimental setup we were able 
to demonstrate that the value of p and the DC action do not de- 
pend on the input state. In fact the AAQPT was realized with 
all the four Bell states entering the DC, obtaining the same 
results of those shown in Fig. [3J}). 

It is worth noting that, even if the AAQPT gives a more 
complete information on the process compared to the imple- 
mented optimal protocol, the latter allows us to achieve a more 
accurate value of p. The inset in FigOa) shows that, for the 
optimal protocol, the measured standard deviation reaches the 
lower bound given by yjp{\ - p)/N [i.e. the square root of Eq. 
(OQ) for the one-dimensional case divided by N], where N is the 
dimension of the sample used to evaluate p, thus demonstrat- 
ing experimentally the attainability of the Cramer-Rao bound. 
We show in the inset in Fig 13b) the standard deviations, well 
above the optimal bound, obtained with the AAQPT (see the 



Supplemental Material for details about the numerical estima- 
tion). The lower optimal bound represented by the black curve 
is below the experimental data, demonstrating that AAQPT is 
far away from the optimal estimation protocol presented in 
this work. 

Conclusion.- An optimal protocol allowing the most ef- 
ficient estimation of a noisy Pauli channel has been experi- 
mentally implemented in this work. The action of the noisy 
channel was introduced on one qubit of a maximally entan- 
gled pair in a controlled way. The efficiency of this method 
has been compared to the one achieved by quantum process 
tomography, demonstrating that the optimal protocol allows 
us to achieve the theoretical lower bound for the errors and to 
perform the estimate of the noisy channel with a lower num- 
ber of measurements. This method can be profitably applied 
when some knowledge on the noise process is available and 
can be successfully implemented in quantum-enhanced tech- 
nologies involving the management of decoherence. 
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SUPPLEMENTARY INFORMATION 

Ancillary Assisted Quantum Process Tomography 
(AAQPT) standard deviations, shown in the inset in Fig. 3b) 
of the main text, have been calculated by a simulation realized 
by MATHEMATICA 5.0. It was based on the following 
procedure: 

• a poissonian uncertainty was associated to the coinci- 
dence counts 



• hundred different matrices of the process were simu- 
lated 

• the fidelity between these matrices and the theoretical 
ones, reported in the main text, was calculated 

• the routine NMAXIMIZE allowed to find numerically 
the value of p maximizing the Fidelity for each simu- 
lated process 



The experimental values of the inset in Fig. 3b) correspond 
to a sample of -1600 coincidences per second. It can be evalu- 
ated that a number of 150000 coincidences per second, which 
is nearly two orders of magnitude larger, is needed to approach 
the standard deviation value obtained by the optimal protocol 
(i.e. the inset in Fig. 3 a) of the main text). 

700 

C/) 600 

1/5 500 
CD 

^ 400- 




FIG. 4. Dip (peak) of the measured coincidence counts as a func- 
tion of the optical path delay for a state \(f)~)(\i//~)) entering the BS in 
absence of noise. 



• this sample, composed of hundred values of p, was used 
to calculate both average and standard deviation. 



